Numerical Analogues of Aronson's Sequence 



Benoit Cloitre 
13 rue Pinaigrier 
Tours 37000, FRANCE 
(Email: abcloitre@wanadoo.fr) 

N. J. A. Sloane 
AT&T Shannon Labs 
Florham Park, NJ 07932-0971, USA 
(Email: njas@research.att.com) 

Matthew J. Vandermast 
53 Piaget Avenue 
Clifton, NJ 07011-1216, USA 
(Email: ghodgesl4@msn.com) 



March 28, 2003 



Abstract 

Aronson's sequence 1, 4, 11, 16, ... is defined by the English sentence "t is the first, 
fourth, eleventh, sixteenth, . . . letter of this sentence." This paper introduces some numerical 
analogues, such as: a(n) is taken to be the smallest positive integer greater than a{n — 1) which 
is consistent with the condition "n is a member of the sequence if and only if a(n) is odd." 
This sequence can also be characterized by its "square", the sequence a^ 2 \n) = a(a(n)), which 
equals 2n + 3 for n > 1. There are many generalizations of this sequence, some of which are 
new, while others throw new light on previously known sequences. 

1. Introduction 

Aronson's sequence, given in the Abstract, is a classic example of a self-referential sequence 
([2], [7], sequence M3406 in [11], A5224 in [10]). It is somewhat unsatisfactory because of the 
ambiguity in the English names for numbers over 100 — some people say "one hundred and 
one", while others say "one hundred one." Another well-known example is Golomb's sequence, 
in which the n th term G(n) (for n > 1) is the number of times n appears in the sequence 
(A1462 in [10]): 

1, 2, 2, 3, 3, 4, 4, 4, 5, 5, 5, 6, 6, 6, 6, 7, 7, 7, 7, 8, . . . 



There is a simple formula for G(n): it is the nearest integer to (and approaches) 




where (j) = (1 + y/h)/2 ([4], [5, Section E25]). 

Additional examples can be found in Hofstadter's books [6], [7] and in [5] and [10]. However, 
the sequence {a(n)} mentioned in the Abstract appears to be new, as do many of the other 
sequences we will discuss. We will also give new properties of some sequences that have been 
studied elsewhere. 

Section 2 discusses the sequence mentioned in the Abstract, and also introduces the "square" 
of a sequence. Some simple generalizations (non-monotonic, "even" and "lying" versions) are 
described in Section 3. The original sequence is based on examination of the sequence mod- 
ulo 2. In Section 4 we consider various "mod y" generalizations. Section 5 extends both the 
original sequence and the "mod y" generalizations by defining the "Aronson transform" of a 
sequence. Finally, Section 6 briefly considers the case when the rule defining the sequence 
depends on more than one term. 

There are in fact a large number of possible generalizations and we shall only mention some 
of them. We have not even analyzed all the sequences that we do mention. In some cases we 
just list the first few terms and invite the reader to investigate them himself. We give the 
identification numbers of these sequences in [10] — the entries there will be updated as more 
information becomes available. 

We have also investigated sequences arising when (2) is replaced by the following rule: 
s(l) = 1, s(n) = s(n — 1) + y if n is already in the sequence, s(n) = s(n — 1) + z otherwise, for 
specified values of x,y,z. This work will be described elsewhere [3]. 

Notation. "Sequence" here usually means an infinite sequence of nonnegative numbers. 
"Monotonically increasing" means that each term is strictly greater than the previous term. 



2. n is in sequence if and only if a(n) is odd 

Let the sequence a(l), a(2), a(3), ... be defined by the rule that a(n) is the smallest positive 
integer > a(n — 1) which is consistent with the condition that 



P = {1,2,3,... 



},N = {0,1,2,3,...}. 



ti, 



n is a member of the sequence if and only if a(n) is odd. 



(i) 
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The first term, a(l), could be 1, since 1 is odd and 1 would be in the sequence. It could 
also be 2, since then 1 would not be in the sequence (because the terms must increase) and 
2 is even. But we must take the smallest possible value, so a(l) = 1. Now a(2) cannot be 2, 
because 2 is even. Nor can a(2) be 3, for then 2 would not be in the sequence but a(2) would 
be odd. However, a(2) = 4 is permissible, so we must take a(2) = 4, and then 2 and 3 are not 
in the sequence. 

So a(3) must be even and > 4, and a(3) = 6 works. Now 4 is in the sequence, so a(4) must 
be odd, and a(4) = 7 works. Continuing in this way we find that the first few terms are as 
follows (this is A79000): 

n: 1 2 3 4 5 6 7 8 9 10 11 12 ■■■ 
a{n): 1 4 6 7 8 9 11 13 15 16 17 18 - 

Once we are past a(2) there are no further complications, a(n — 1) is greater than n, and 

we can, and therefore must, take 

a(n) = a(n- 1) +e , (2) 

where e is 1 or 2 and is given by: 

a(n — 1) even a(n — 1) odd 

n in sequence 1 2 

n not in sequence 2 1 

The gap between successive terms for n > 3 is either 1 or 2. 

The analogy with Aronson's sequence is clear. Just as Aronson's sentence indicates exactly 
which of its terms are t's, {a(n)} indicates exactly which of its terms are odd. 

We proceed to analyze the behavior of this sequence. 

First, all odd numbers > 7 occur. For suppose 2t + 1 is missing. Therefore a{i) = 2t, 
a{i + 1) = 2t + 2 for some i > 3. From the definition, this means i and i + 1 are missing, 
implying a gap of at least 3, a contradiction. 

Table I shows the first 72 terms, with the even numbers colored red. 

Examining the table, we see that there are three consecutive numbers, 6, 7, 8, which are 
necessarily followed by three consecutive odd numbers, a(6) = 9, a(7) = 11, a(8) = 13. Thus 9 
is present, 10 is missing, 11 is present, 12 is missing, and 13 is present. Therefore the sequence 
continues with o(9) = 15 (odd), a(10) = 16 (even), a(13) = 19 (odd), a(14) = 20 (even). 
This behavior is repeated for ever. A run of consecutive numbers is immediately followed by 
a run of the same length of consecutive odd numbers. 
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Table I: The first 72 terms of the sequence "n is in sequence if and only if a(n) is odd. 



n : 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


a(n) : 


1 


4 


6 


7 


8 


9 


11 


13 


15 


16 


n : 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


a{n) : 


17 


18 


19 


20 


21 


23 


25 


27 


29 


31 


n : 


21 


22 


23 


24 


25 


26 


2i 


28 


29 


30 


a(n) : 


33 


34 


35 


36 


37 


38 


39 


40 


41 


42 


n : 


31 


32 


33 


34 


35 


36 


37 


38 


39 


40 


a(n) : 


43 


44 


45 


47 


49 


51 


53 


55 


57 


59 


n : 


41 


42 


43 


44 


45 


46 


47 


48 


49 


50 


a{n) : 


61 


63 


65 


67 


69 


70 


71 


72 


73 


74 


n : 


51 


52 


53 


54 


55 


56 


57 


58 


59 


60 


a(n) : 


75 


76 


77 


78 


79 


80 


81 


82 


83 


84 


n : 


61 


62 


63 


64 


65 


66 


67 


68 


69 


70 


a(n) : 


85 


86 


87 


88 


89 


90 


91 


92 


93 


95 



n : 71 72 • • • 

a(n) : 97 99 

Let us define the /c th segment (for k > 0) to consist of the terms a(n) with n = 9 • 2 k — 3 + j 
where —3 • 2 k < j < 3 • 2 k — 1. In the table the segments are separated by vertical lines. The 
first half of each segment, the terms where j < 0, consists of consecutive numbers given by 
a(n) = 12 • 3 fc — 3 + j; the second half, where j > 0, consists of consecutive odd numbers given 
by a(n) = 12 • 3 fc - 3 + 2j. We can combine these formulae, obtaining an explicit description 
for the sequence: 

o(l) = l, o(2)=4, 

and subsequent terms are given by 

a(9 • 2 k - 3 + j) = 12 • 2 k - 3 + |? + (3) 

for k > 0, -3- 2 k < j < 3- 2 k . 

The structure of this sequence is further revealed by examining the sequence of first differ- 
ences, Aa(n) = a(n + 1) — a{n), n > 1, which is 

3,2,1,1,1,2,2,2,1 6 , 2 6 ,1 12 , 2 12 ,1 24 ,2 24 ,... (4) 
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(A79948), where we have written l m to indicate a string of m l's, etc. The oscillations double 
in length at each step. 

Segment begins with an even number, 6, but all other segments begin with an odd 
number, 9 • 2 k — 3. All odd numbers occur in the sequence except 3 and 5. The even numbers 
that occur are 4, 6, 8 and all numbers 2m with 

9 . 2 fe - 1 - 1 < m < 6 • 2 fe - 2, k > 1 . 

The sequence of differences, (4), can be constructed from the words in a certain formal 
language (cf. [9]). The alphabet is {1,2,3} and we define morphisms 0(1) = 2,9(2) = 1,1. 
Then (4) is the concatenation 

S-i, Sq, S\, S2, ■ ■ ■ , (5) 

where 

5_i = {3, 2}, S = {1, 1, 1}, S k+ i = 9(S k ) for k > . (6) 

To prove this, note that for n > 3, a difference of 2 only occurs in {a(n)} between a pair 
of odd numbers. Suppose a(i) = 2j + 1, a(i + 1) = 2j + 3; then a(2j + 1) = 2x + 1 (say), 
a(2j + 2) = 2x + 2, a(2j + 3) = 2x + 3, producing two differences of 1. Similarly, if there is a 
difference of 1, say a(i) = j, a(i + 1) = j + 1, then a(j) = 2x + 1, a(j + 1) = 2x + 3, a difference 
of 2. 

The ratio n/a(n), which is the fraction of numbers that are in the sequence, rises from close 
to 2/3 at the beginning of segment k (assuming k is large), reaches a maximum 3/4 at the 
midpoint of the segment, then falls back to 2/3 at the end of the segment. It is not difficult to 
show that if n is chosen at random in the k th segment then the average value of the fraction 
of numbers in the sequence at that point approaches 

3_i log ^ = 0.7075... 
4 4 & 27 

for large k. 

The sequence has an alternative characterization in terms of its "square." 
The square of a sequence s = {s(n) : n > uq} is given by s^ 2 ) = {s(s(n)) : n > no}. If s is 
monotonically increasing so is . 

Lemma 1. Let s be monotonically increasing. Then n (> no) is in the sequence s if and only 
if s(n) is in the sequence s^ 2 \ 
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Proof. If n is in the sequence, n = s(i) for some i > no, and s(n) = s(s(i)) is in s^ 2 ). 
Conversely, if s(n) G s^ 2 \ s(n) = s(s(i)) for some i > no, and since s is monotonically 
increasing, n = s(i). ■ 

For our sequence a = {a(n)}, examination of Table I shows that a^ 2 ) = {1, 5, 7, 9, 11, . . .} = 
{1} U 2P + 3. This can be used to characterize a. More precisely, the sequence can be defined 
by: a(l) = 1, a(2) = 4, a(3) = 6 and, for n > 4, a(n) is the smallest positive integer which is 
consistent with the sequence being monotonically increasing and satisfying a(a(n)) = 2n + 3 
for n > 2. 

This is easily checked. Once the first three terms are specified, the rule a(a(n)) = 2n + 3 
determines the remaining terms uniquely. 

In fact that rule also forces a(2) to be 4, but it does not determine a(3), since there is an 
earlier sequence {a'(n)} (in the lexicographic sense) satisfying a'(l) = 1, a'(a'(n)) = 2n + 3 for 
n > 2, namely 

1,4, 5, 7, 9, 10, 11, 12, 13, 15, 17, 19, 21, 22, ... , 
(A80596), and given by a'(l) = 1, 

«'(6 - 2 fc - 3 + j) = 8 • 2 fc - 3 + |i + (7) 

for k > 0, -2 fe+1 < j < 2 k+1 . 

As the above examples show, the square of a sequence does not in general determine the 
sequence uniquely. A better way to do this is provided by the "inverse Aronson transform", 
discussed in Section 5. 

3. First generalizations 

The properties of {a(n)} given in Section 2 suggest many generalizations, some of which will 
be discussed in this and the following sections. 

(3.1) Non-monotonic version. If we replace "a(n) > o(n — 1)" in the definition by "a(n) 
is not already in the sequence", we obtain a completely different sequence, suggested by J. C. 
Lagarias (personal communication): 6(n), n > 1, is the smallest positive integer not already 
in the sequence which is consistent with the condition that "n is a member of the sequence if 
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and only if b(n) is odd." This sequence (A79313) begins: 

n: 123456789 10 
b(n) : 1 3 5 2 7 8 9 11 13 12 

n : 11 12 13 14 15 16 17 18 19 20 
b(n) : 15 17 19 16 21 23 25 20 27 29 

The even members are shown in red. The behavior is simpler than that of {a(n)}, and we 
leave it to the reader to show that, for n > 5, b(n) is given by 

b(At-2) = At, 

b(At - 1) = 6t - 3 , 

b{At) = 6t-l, 

6(4t + l) = 6t + l. 

All odd numbers occur. The only even numbers are 2 and At, t > 2. (The square 

b (2) 

is not 

so interesting.) 

(3.2) "Even" version. If instead we change "odd" in the definition of {o(n)} to "even", we 
obtain a sequence c which is best started at n = 0: c(n), n > 0, is the smallest nonnegative 
integer > c(n — 1) which is consistent with the condition that 

"n is a member of the sequence if and only if c(n) is even." (8) 

This is A79253: 0, 3, 5, 6, 7, 8, 10, 12, 14, 15, .... It is easily seen that c(n) = a(n + 1) -1 for 
n > 0, so there is nothing essentially new here. Also c^ 2 ) = {0} U 2P + 4. 

(3.3) The "lying" version. The lying version of Aronson's sequence is based on the com- 
pletely false sentence "t is the second, third, fifth, sixth, seventh, . . . letter of this sentence." 
The sentence specifies exactly those letters that are not t's, and produces the sequence (A81023) 
2, 3, 5, 6, 7, 8, 9, 10, 11, 12, .... 

Just as {a(n)} is an analogue of the original sequence, we can define an analogue {d(n) : 
n > 1} of this sequence by saying that: d(n) is the smallest positive integer > d(n — 1) such 
that the condition "n is in the sequence if and only if d(n) is odd" is false. Equivalently, the 
condition "either n is in the sequence and d(n) is even or n is not in the sequence and d(n) is 
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odd" should be true. The resulting sequence (A80653) begins 2, 4, 5, 6, 8, 10, 11, 12, 13, 14, 
.... We will give an explicit formula for d(n) in the next section. 

A related sequence is also of interest. Let {d'(n)} be defined by d'(l) = 2, and, for n > 1, 
d'(n) is the smallest integer greater than d'(n — 1) such that the condition "n and d'(d'(n)) 
have opposite parities" can always be satisfied. One can show that this is the sequence 

2,4,5,7,8,9,11,12,13,14, 16,... 

(A14132), the complement of the triangular numbers, with d'(n) = n + nearest integer to 

4. The "mod m" versions 

Both {a(n)} and {c(n)} are defined modulo 2. Another family of generalizations is based on 
replacing 2 by some fixed integer y > 2. To this end we define a sequence {s(n) : n > no} by 
specifying a starting value s(rao) = so> an d the condition that n is in the sequence if and only 
if s(n) = z (mod y), where y and z are given. 

Although we will not digress to consider this here, it is also of interest to see what happens 
when "if and only if" in the definition is replaced by either "if" or "only if." (We mention 
just one example. The above sequence {d(n)}, prefixed by d(0) = 0, can be defined as follows: 
d(n) is the smallest nonnegative number > d(n — 1) such that the condition "n (n > 0) is in 
the sequence only if d(n) is even" is satisfied.) 

We saw in the previous section that {a(n)} can also be characterized by the property that 
its square a^ 2 \n) = a(a(n)) is equal to 2n + 3 for n > 2 (together with some appropriate 
initial conditions). This too can be generalized by specifying that the sequence {s(n)} satisfy 
s(s(n)) = yn + z, for given values of y and z. The two generalizations are related, but usually 
lead to different sequences. The s(s(n)) family of generalizations will connect the present 
investigation with several sequences that have already appeared in the literature. There are 
too many possibilities for us to give a complete catalogue of all the sequences that can be 
obtained from these generalizations. Instead we will give a few key examples and one general 
theorem. Many other examples can be found in [10]. 

A simple "mod 3" generalization is: e(l) = 2, and, for n > 1, e(n) is the smallest integer 
> e(n — 1) which is consistent with the condition that 

"n is a member of the sequence if and only if e(n) is a multiple of 3." (9) 
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This turns out to be James Propp's sequence 

2,3,6,7,8,9,12,15,18,19,... , 

which appeared as sequence M0747 in [11] (A3605in [10]). Propp gave a different (although 
equivalent) definition involving the square of the sequence: {e(n)} is the unique monoton- 
ically increasing sequence satisfying e(e(n)) = 3n for all n > 1. Michael Somos (personal 
communication) observed that this sequence satisfies 

e(3n) = 3e(n), 
e(3n + l) = 2e(n) + e(n + 1), 

e(3n + 2) = e(n) + 2e(n + 1). (10) 

An analysis similar to that for {a(n)} leads to the following explicit formula, which appears to 
be new: 

e(2-3 fc + i) = 3 fc+1 + 2j + |j|, (11) 

for k > and -3 fc < j < 3 k . 

A sequence closely related to {e(n)} had earlier been studied by Arkin et al. [1, Eq. (12)]. 
This is the sequence e'(n) = e(n) — n (A6166). Arkin et al. give a recurrence similar to (10). 

The sequence {e(n)} can be generalized as follows. 

Theorem 1. Let y and z be integers of opposite parity satisfying 

y>2, y + z> 1, 2y + z>4. (12) 

Then there is a unique monotonically increasing sequence {f(n)} satisfying /(l) = ^(y + z+1) 
and f(f(n)) = yn + z for n > 1. It is given by 

.k * , A _ ,.k+i z V + !, , y- !, 



/or fc > 0, where 



and 



y + z — 1 k ^ . y + z—1 k 

- — y k < j < - — y k 



(y + l)(y + z-l) y + z-1 

2(2/ — i) y-i 
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Proof. The sequence begins with /(l) = (y + z + l)/2, and is constrained by f(f(l)) = 
/((y + z + l)/2) = y + z. Since it is monotonically increasing, the intermediate terms are 
forced and are given by f(l + i) = (y + z + l)/2 + i for 0<i<(y + z + l)/2. The terms 
/(2), /(3), . . . , /((y + z + l)/2) now determine the values of /(/(2)) = /(y + z), /(/(3)) = 
/(2y + z), . . ., and we find that the stretch from /((y + z + l)/2) through /(y + z) is given by 
/((y -+- z + l)/2 + i) = y + z + yi for < i < (y + z — l)/2. Continuing in this way we find 
that the sequence is completely determined, and (13) follows after relabeling the indices. The 
conditions (12) ensure that there is no contradiction in calculating the initial terms, and once 
started the sequence has a unique continuation. ■ 

It can be shown (we omit the details) that this sequence can also be defined by: f(l) = 
(y + z + l)/2, and, for n > 1, f(n) is the smallest integer > f(n — 1) which is consistent with 
the condition that "n is a member of the sequence if and only if /(n) belongs to the set 

[2,...,i(y + z-l)] U {iy + z:i > 1} ." (14) 

If (y + z — l)/2 < 1, the first set in (14) is to be omitted. 

Examples. Setting y = 3, z = in the theorem produces {e(n)}. 

Setting y = 2, z = 1 yields another interesting "mod 2" sequence. This is the sequence 
{g(n) : n > 1} that begins 

2,3,5,6,7,9,11,12,13, 14,... 

(A80637). It has the following properties: 

(i) By definition, this is the unique monotonically increasing sequence {g(n)} satisfying 
g(l) = 2, g(g(n)) = 2n + 1 for n > 2. 

(ii) n is in the sequence if and only if g(n) is an odd number > 3. 

(iii) The sequence of first differences is (A79882): 

1,2,1 2 ,2 2 ,1 4 ,2 4 ,1 8 ,2 8 ,1 16 ,2 16 ,.... 

(iv) 

5 (3 • 2 k - 1 + j) = 2 • 2 k+l - 1 + |? + , 

for fc > 0, -2 k <j<2 k (from (13)). 

(v) g(2n) = g(n) + y(n - 1) + 1, y(2n + 1) = 2g(n) + 1, for n > 1 (taking y(0) = 0). 
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(vi) The original sequence {a(n)} satisfies a(3n) = 3g(n), a(3n+l) = 2g(n)+g(n+l), a(3n+ 
2) = g( n ) + 2#(n + 1), for n > 1. 

(vii) The "lying version" of Section 3 is given by d{n) = g(n + 1) — 1 for n > 1. 

(viii) Let </(n) = g(n) + 1. The sequence {g'(n) : n > 2} was apparently first discovered by 
C. L. Mallows, and is sequence M2317 in [11] (A7378 in [10]). This is the unique monotonically 
increasing sequence satisfying g'(g'{n)) = 2n. An alternative description is: g'(n) (for n > 2) 
is the smallest positive integer > g'{n — 1) which is consistent with the condition that 

"n is a member of the sequence if and only if g'(n) is an even number > 4". (15) 

Note that, although (8) and (15) are similar, the resulting sequences {c(n)} and {g'(n)} are 
quite different, g' is not directly covered by Theorem 1, and we admit that we have not been 
able to identify the largest family of sequences which can be described by formulae like (3), 
(7), (11), (13). 

The sequence {h(n)} defined by: /i(l) = 2, and, for n > 1, h{n) is the smallest positive 
integer > h(n — 1) which is consistent with the condition that "n is a member of the sequence 
if and only if h{n) is a multiple of 6": 

2,6,7,8,9,12,18,24,30,31,... , 

(A80780), shows that such simple rules do not hold in general. We can characterize the 
sequence of first differences in a manner similar to (5), (6): the alphabet is {1, 2, . . . ,6}, and 
we define morphisms 9(i) = 1, 1, . . . , 1, 7 — i (with i — 1 l's followed by 7 — i), for i = 1, . . . , 6. 
Then the sequence of differences of {h(h)} is So, Si, S2, . . ., where So = {4}, S^+i = 6(Sk) for 
k > 0. However, it appears that no formula similar to (3) holds for h{n). 

We end this "mod m" section with two interesting "mod 4" sequences. The even numbers 
satisfy s(s(n)) = 4n, and the odd numbers satisfy s(s(n)) = 4n + 3. But there are lexicograph- 
ically earlier sequences with the same properties. The "fake even numbers" {i(n) : n > 0} are 
defined by the property that i(n) is the smallest nonnegative integer > i(n — 1) and satisfying 
i(i(n)) = An (A80588): 

0,2,4,5,8, 12,13,14, 16,17,... 
We analyze this sequence by describing the sequence of first differences, which are 

2,2,1,3,4,1,1,2,1,1,1,1,4,4,1,3,... 
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After the initial 2, 2, 1, this breaks up into segments of the form 



3 Sk 2 T k 



where T/% is the reversal of 



l 1 4 2 l 4 4 8 l 16 4 32 



...4 




and Sk is the reversal of 



j2 4 1 X 8 4 2 j_32 ^e...^" 1 ^ 



,2fc— 2 



The "fake odd numbers", i'(n), are similarly defined by i'(i'(n)) = 4n + 3: 1, 3, 4, 7, 11, 12, 
13, 15, 16, 17, ... (A80591), and satisfy i'(n) = i(n + 1) - 1. 

5. The Aronson transform 

A far-reaching generalization of both the original sequence and the "mod m" extensions of the 
previous section is obtained if we replace "odd number" in the definition of {a(n)} by "member 
of /3" , where (3 is some fixed sequence. 

More precisely, let us fix a starting point no, which will normally be or 1. Let (3 = {/3(n) : 
n > n o} be an infinite monotonically increasing sequence of integers > no with the property 
that its complement (the numbers > no that are not in (3) is also infinite. Then the sequence 
a = {a(n) : n > no} given by: a(n) is the smallest positive integer > a(n — 1) which is 
consistent with the condition that 



is called the Aronson transform of (3. 

Theorem 2. The Aronson transform exists and is unique. 

Proof. For ease of discussion let us call the numbers in (3 "hot" , and those in its complement 
"cold." We will specify the transform a, leaving to the reader the easy verification that this 
has the desired properties, in particular that there are no contradictions. 

The proof is by induction. First we consider the initial term a(no). If no is hot, a(no) = no- 
If no is cold, a(no) = smallest cold number > no + 1. 

For the induction step, suppose a(n) = k for n > hq. 



LL. 



; n is in a if and only if a(n) is in (3' 
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Case (i), k = n. If n + 1 is hot then a{n + 1) = n + 1. If n + 1 is cold then a(n + 1) = 
smallest cold number > n + 2. 

Case (ii), k > n. If /c = n + 1 then a(n + 1) = smallest hot number > n + 2. If 
/c > n + 1 then if + 1 is hot, a(n + 1) = smallest hot number > k + 1, while if re + 1 is cold, 
a{n + 1) = smallest cold number > k + 1. ■ 

In certain cases it may be appropriate to specify some initial terms in a to get it started 
properly. 

Examples. Of course taking j3 to be the odd numbers (with hq = 1) leads to our original 
sequence {a(n)}, and the even numbers (with no = 0) lead to {c(n)} of Section 3. 

The sequences P (with hq = 1) and N (with no = 0) are fixed under the transformation. 

If we take /3 to be the triangular numbers we get 1, 4, 5, 6, 10, 15, 16, 17, 18, 21, . . . 
(A79257); the squares give 1, 3, 4, 9, 10, 11, 12, 13, 16, 25, . . . (A79258); the primes give 4, 6, 
8, 11, 12, 13, 14, 17, 18, 20, ... (A79254); and the lower Wythoff sequence (A201), in which 
the n th term is [n<j)\, gives 1, 5, 7, 10, 11, 13, 14, 15, 18, 19, . . . (A80760). 

Taking the Aronson transform of {a(n)} itself we get 1, 3, 4, 6, 10, 11, 12, 14, 22, 23, . . . 
(A79325). [10] contains several other examples. 

The inverse transform may be defined in a similar way. Given an infinite monotonically 
increasing sequence a. = {a(n) : n > no} of numbers > no, such that its complement (the 
numbers > no that are not in a) is also infinite, its inverse Aronson transform is the sequence 
P = {f3(n) : n > hq} such that the Aronson transform of (3 is ex. 

Theorem 3. The inverse Aronson transform exists and is unique. 

Proof. We establish this by giving a simple algorithm to construct the inverse transform. We 
illustrate the algorithm in Table II by applying it to the sequence of squares, cx = {n 2 : n > 0}. 

Form a table with four rows. In the first row place the numbers n = no, no + 1, no + 2, . . ., 
and in the second row place the sequence a (no), a (no + 1), . . .. The third row contains what 
we will call the "hot" numbers: these will comprise the elements of the inverse transform. The 
fourth row are the "cold" numbers, which are the complement of the hot numbers. 

The third and fourth rows are filled in as follows. If n is in (resp. not in) the sequence a, 
place a{n) in the n-th slot of the hot (resp. cold) row. 
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To complete the table we must fill in the empty slots. Suppose we are at column n, where 
we have placed a(n) in one of the two slots. Let l n be the largest number mentioned in columns 
no, — 1 in the hot or cold rows. Then we place the numbers l n + 1, . . . , a(n) — 1 in the 
empty slot in column n, with the single exception that if n is not in the sequence and l n = n— 1 
then we place n in the cold slot rather than the hot slot. (This is illustrated by the position 
of 2 in the fourth row of Table II.) We leave it to the reader to verify that the entries in the 
"hot" row form the inverse Aronson transform (3. m 

Table II: Computation of inverse Aronson transform of the squares. The "hot" numbers com- 
prise the transform. 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


a n 


1 


4 


9 


16 


25 


26 


49 


64 


81 


"hot" 


1 


3 


5-8 


16 


17-24 


26-35 


37-48 


50-63 


81 


"cold" 




2,4 


9 


10-15 


25 


36 


49 


64 


65-80 



It follows from Lemma 1 that (3 contains the members of a^ 2 \ but in general (3 ^ a^ 2 \ 
The additional terms in (3 serve to make it possible to recover a uniquely from (3. 

Examples. As shown in Table II, the inverse Aronson transform of the squares (A10906) is 

1; 3; 5, 6, 7, 8; 16, 17,... , 24; 26, 27, . . . . 

This consists of a number of segments (separated here by semicolons). For k > 1 the /c th 
segment is {k 2 } if A; is a square, or {(k — l) 2 + 1, . . . , k 2 — 1} if k is not a square, except that 
the second segment = {3}. 

The inverse transform of the primes is 3, 5, 6, 11, 12, 17, 18, 20, 21, 22, ... (A80759) - 
this has a similar decomposition into segments. 

The inverse transform of the lower Wythoff sequence is 1, 4, 6, 7, 9, 10, 12, 14, 15, 17, . . . 
(A80746). This consists of the numbers [<f>k\ +k — l(k>l) and [2#J + k - 1 (k > 2). The 
inverse transform of our original sequence {a(n)} is the sequence of odd numbers (whereas, as 
we saw in Section 2, omits 3). 

In general (because of the above algorithm), the inverse Aronson transforms are easier to 
describe than the direct transforms. 
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6. More complicated conditions 

Finally, we may make the condition for n to be in the sequence depend on the values of several 
consecutive terms a(n), a(n + 1), . . . , a(n + r), for some fixed r. To pursue this further would 
take us into the realm of one-dimensional cellular automata (cf. [8], [12]), and we will mention 
just two examples. 

q(n) is the smallest positive integer > q(n — 1) which is consistent with the condition that 
"n is in the sequence if and only if q{n) is odd and q(n — 1) is even" (A79255): 

1,4,6,9,12,15,18,20,23,26,28,... 

The gaps between successive terms are always 2 or 3. Changing the condition to ". . . both 
q{n) and q(n + 1) are odd" gives A79259: 

1,5,6,10,11,15,19,20,24,25,... 
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